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1 .  Introduction 

The  background  for  most  of  the  work  described  in  this  report  has  been 
given  previously  and  no  attempt  is  made  to  repeat  it  here.  The  reader  who 
wants  further  details  should  refer  to  the  original  or  continuation  proposals 
entitled  "Correlation  and  Collective  Modes  in  Narrow  Band  Materials:  NiO 
and  FeO,"  under  Grant  No.  AFOSR-80-0023 .  The  progress  in  this  research 
during  the  past  year  is  described  in  Section  III.  An  important  new  component 
of  our  studies  of  electronic  excitations  and  collective  modes  is  work  on  the 
effect  of  plasmons  and  excitons  on  superconductivity  in  transition  metal 
oxides  which  is  described  in  Section  1I1-C. 
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I  I .  Research  Objectives 

The  overall  objective  of  this  research  is  to  study  the  effects  of 
electronic  excitations,  especially  plasmons  and  interband  excitations,  on  the 
properties  of  materials  for  which  these  effects  are  significant.  The  materials 
we  have  chosen  are  narrow  band  systems,  such  as  the  transition  metal  oxides. 

For  these  materials,  it  is  necessary  to  emphasize  their  localized  properties 
and,  especially,  strong  electron-lattice  interactions  in  constructing  the 
correlation  and  response  functions.  In  the  first  year  of  this  work,  we  have 
set  up  this  description  of  the  electronic  structure  in  terms  of  a  local 
orbitals  approach  and  have  begun  some  calculations. 

Our  studies  of  electronic  correlation  and  its  effect  on  the  band 
structure  have  begun  with  TiO.  We  will  continue  with  investigations  of  systems 
where  the  effects  of  electronic  excitations  are  even  stronger,  such  as  FeO  and 
NiO.  For  these  systems  interband  excitations  become  quite  prominent  and 
complicate  the  optical  or  electron  energy  loss  spectra.  We  will  also  continue 
our  studies  of  collective  modes  for  these  systems,  such  as  plasmons,  in  the 
presence  ot  stronu  lattice  interactions. 

We  are  investigating  the  contributions  of  electronic  excitations  to 

quantities  which  can  be  measured  experimentally.  Presently  we  are 

analyzing  the  electronic  loss  structure  of  TiO  and  FeO.  In  the  future  we 

will  extend  this  analysis  to  NiO  and  other  narrow  band  systems.  We  will 

want  to  know  how  the  loss  structures  due  to  plasmon  and  interband 

excitations  behave  in  these  materials,  and  we  will  want  to  develop  methods  for 

explaining  the  significant  changes  that  occur  in  some  cases  when  the  surface 

is  contaminated  or  the  material  is  implanted  with  impurities.  From  these 

studies  we  hope  to  develop  an  important  tool  for  explaining  the  effects  of 

surlace  or  bulk  contamination  on  these  systems. 

AIR  FORCE  OFFITO  OF  SCIENTIFIC  RESEARCH  (AISC) 
NOTICE  OF  TRANSMITTAL  TO  POC 

This  technical  report  has  bean  reviewed  and  is 
approved  for  public  release  IAW  AFR  190-12  (7b). 
Distribution  is  unlimited. 

A.  D.  BLOSF 

_ _ _  Technical  Information  Officer 
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We  are  investigating  plasmon  and  exciton  coupling  in  the  establishment  of 
superconductivity  for  narrow  band  systems.  We  will  want  to  know  how  electrons 
coupled  by  these  mechanisms  shield  external  fields,  and  what  effect  the  strong 
coupling  has  on  the  superconducting  properties.  We  have  begun  studies  on 
Cu„0  lor  which  some  interesting  experimental  evidence  exists.  In  the  future  we 
will  move  to  other  transition  metal  oxides. 
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r II.  Progress 

A.  Correlation  in  Narrow  Band  Materials. 

The  effect  of  correlation  on  the  band  structure  of  a  material  is  contained 
in  the  self-energy  matrix  £.  To  find  the  correlated  energy  bands  and  wave 
functions,  one  solves  the  energy  dependent  matrix  equation 


So  +  I<“) 


-*  -► 
A(w)  =  E(co)  A(w) 


(3.1) 


where  Hq  is,  for  example,  a  one-  actron  Hamiltonian  such  as  Hartree  Fock  (HF) 
and  u>  is  the  frequency.  Now  in  general  we  must  diagonalize  Eq.  (1)  to  find 
the  energy  bands  E(w).  However,  if  the  diagonal  terms  in  Hq  +  £  are  dominant, 
the  correlated  energy  bands  are  just 


k  -  e!f  +  iV.ej  »  (3-2> 

k  k  k 

assuming  that  Hq  is  the  HF  Hamiltonian.  The  self-energy  to  second  order  in  the 
coulomb  interaction  between  Bloch  states  is 
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co  +  E  :*•,  -  E  —  E  „  • ,  “► 

mk  Jlk'+q  l  k+q 


where  (m'  ►  ?')  signifies  that  the  second  set  of  V's  is  the  same  as  the  first 
except  that  m'  is  replaced  by  l ' .  Here  V  refers  to  the  coulomb  interaction, 
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and  m  always  refers  to  particle  Bloch  states  (E  >  E„)  while  £  always  refers 

m  r 

to  hole  states  (E^,  <  E^). 

There  is  a  great  savings  in  computational  effort  to  be  had  by  converting 
Eq.  (3.3)  into  a  sum  over  Wannier  states  (for  details  see  proposal  for 
continuation  of  AFOSR-8 0-0023) .  This  savings  makes  the  solutions  of  Eqs.  (3.1- 
3.2)  possible.  The  transformation  of  Eq.  (3.3)  to  the  Wannier  representation 
requires  use  of  the  relation 


V'-V 


7=  l 

✓N  BZ 


,-k  ik-R 
*nk(r)e  1 


(3.4) 


between  the  Wannier  and  Bloch  states  where  N  is  the  number  of  unit  cells  and 
n  is  the  band  index.  Thus  we  must  remove  factors  like  e  j  from  the  coulomb 
matrix  elements  in  Eq.  (3.3)  so  that  we  have  the  sum 
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(3.4*) 


Then  the  V's  are  matrix  elements  over  the  Wannier  states.  We  see  that  the  sums 
over  band  index  in  Eq.  (3.3)  are  limited  to  states  above  or  below  the  Fermi 
level,  and  this  restriction  will  reduce  computational  effort.  More  importantly 
our  Wannier  functions  are  restricted  to  the  band,  complex  which  contributes  most 
strongly  to  correlation  in  a  given  energy  region  of  the  band  structure.  We  will 
see  that  this  group  of  bands  can  be  much  smaller  than  the  total  number  of  bands. 

The  effect  of  correlation  on  the  valence  bands  of  NiO  and  FeO  is  still  not 
fully  explained.  Especially  for  NiO,  the  detailed  structure  and  width  of  the 
d-band  region  is  uncertain.  We  are  planning  on  doing  calculations  of  correlation 
corrections  for  these  two  materials,  but  first  we  have  undertaken  a  calculation 
of  72  in  Eq.  (3.3)  for  TiO.  This  material  was  the  subject  of  some  controversy 
a  few  years  ago.  An  ab-lnltlo  Hartree  Fock  calculation  by  Jennison  and  Kunz^ 
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sliowed  hybridized  oxygen  2p  and  titanium  3d  bands  with  no  energy  gap  between 

2 

them.  These  results  differed  sharply  from  local  exchange  calculations  which 

show  a  band  gap  of  at  least  2eV  between  the  oxygen  2p  and  titanium  3d  bands. 

3 

Experimental  X-ray  photoemission  results  [by  Ichikawa  et  al.  ]  for  what  was 

evidently  TiO  seemed  to  support  the  Jennison  and  Kunz  calculation.  However, 

4 

later  experiments  showed  that  the  Ichikawa  samples  must  have  oxidized  to 

become  TiO^  and  thus  the  d  band  structure  in  the  XPS  spectrum  was  missing. 

4 

The  correct  TiO  spectra  had  a  band  gap  and  were  in  agreement  with  the 
local  density  calculations.  Subsequently,  Kunz^  argued  that  the  gap  was  due 
to  correlation  corrections  neglected  in  the  Hartree  Fock  calculation  and 
asserted  that  relaxation  would  raise  the  d  bands  with  respect  to  the  p  bands. 

Our  calculation  should  test  Kunz'  assertion  and  will  also  provide  an  assessment 
of  correlation  effects. 

We  have  included  9  bands  in  the  TiO  self-energy  calculation  which  arise 
from  the  titanium  4s  and  3d  states  and  the  oxygen  2p  states.  The  other 
bands  are  well  out  of  the  valence  -  conduction  band  region,  and  we  consider  them 
to  not  contribute  much  to  the  self-energy  near  the  Fermi  surface.  Except  for  a 
small  region  around  the  T  point  in  the  uncorrelated  band  structure  where  all  the 
bands  are  above  the  Fermi  level,  three  bands  are  occupied,  three  are  unoccupied, 
and  three  are  partly  occupied.  Thus  by  working  in  the  Wannier  representation, 
we  have  an  immediate  partioning  and  need  only  consider  6  bands  for  the  occupied 
states  and  6  bands  for  the  unoccupied  states. 

The  self-energy  calculation  still  requires  a  major  computational  effort 
in  spite  of  the  partitioning.  We  must  evaluate  the  self-energy  matrix  at  k-points 
throughout  the  Brillouin  zone  and  at  frequencies  w  which  cover  the  valence  band 
energies  for  each  k-point.  In  the  TiO  calculation  we  have  chosen  10  frequency 
points  at  each  k  to  evaluate  £(k,w)  and  must  interpolate  between  the  solutions 
of  Kq.  (3.1)  to  find  En=E(En)  for  a  given  band  n.  We  have  used  every  possible 
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symmetry  to  make  the  calculation  manageable.  The  computation  of  the  Wannier 
coef  f icients 


C 

na 


->• 

ik*  Ri 


(3.5) 


is  greatly  reduced  by  the  fact  that  all  rotations  of  the  lattice  set  are 

equivalent  and  thus  we  need  calculate  only  one  (^(R^)  for  each  shell  of 

atoms.  Futhermore,  since  our  system  has  inversion  symmetry,  we  can  choose 
-k  k* 

C  -  C  and  thus  Eq.  (3.5)  becomes 
na  na  ^ 


C 

na 


(Rl)  =  l  Re 
k 


Ck  eik‘M 
na  I 


(3.6) 


The  index  a  normally  runs  over  all  the  basis  functions  in  the  basis  set 

(e.g.,  40  functions  for  TiO).  This  number  of  functions  can  be  reduced  to  9, 

the  number  of  bands  in  our  complex,  by  variationally  choosing  optimized  basis 

functions.3  However,  since  in  this  case  we  are  using  local  orbitals  which 

are  in  a  sense  already  optimized  for  each  band,  we  may  simply  choose  the  9  local 

orbitals  contributing  most  to  the  Bloch  functions  in  our  band  complex.  Once 

calculated,  the  C  (Ri)  is  used  to  construct  the  coulomb  matrix  elements 
na  1 

V  „  _  which  appear  with  expression  (3.4')  in  the  construction 

n,Rl,n2R2;  n^.n^ 

of  y2(k,w).  If  we  suppress  the  band  indices  n  and  basis  functions  indices  a, 
we  may  write  the  coulomb  matrix  element  as 


V  if  o  =  l  C(R  -R  )  C(R’-R  )  C(R’-R-) 

W  J  °  R, R_R  R.  11  22  33 

12  3  4 


C(5,,-R.)  V-+  -+  ►  -V  ►  ± 

4  4^  R.-R, ,R„-R. ;  R.-R. ,( 

1  4  2  4  3  4 


(3.7) 


x 


showing  that  we  need  only  consider  interactions  which  include  the  origin  in 

our  sum.  Since  the  sum  in  Eq.  (3.3)  involves  only  certain  combinations  of 

occupied  and  unoccupied  bands,  not  all  matrix  elements  V'  need  be  evaluated. 

v2 

The  final  term  appearing  in  )  is  the  expression  appearing  in  (3.4'). 
This  expression  is  similar  in  structure  to  the  free  electron  polarization 
and  thus  we  call  it  x(k;R^,R2)*  We  may  write  (3.4’)  as 


x(k;R1,R2) 


l  x(k,k ',q ,;R1,R2) 


(3.8) 


again  suppressing  all  band  indices.  To  calculate  every  element  in  Eq.  (3.8) 
would  require  a  large  amount  of  time  and  a  massive  amount  of  storage. 
Fortunately,  not  every  element  need  be  calculated  since  many  are  related  by 
symmetry.  The  basic  relation  is 


X(Rk,Rk\Rq  ;  Rltit2)  =  x(k,k\q  ;  R-1^  R  ^  (3.9) 

where  R  is  a  rotation  belonging  to  the  rotation  group  of  the  crystal.  In 

particular  Eq.  (3.9)  indicates  that  we  need  calculate  x  for  only  those  k 

values  in  an  irreducible  wedge  of  the  Brillouin  zone  if  we  include  all 

rotations  of  the  other  variables.  Furthermore  all  rotation  of  the  type 
.  •+  — ►  ■  ► 

R  k  =  Rk  where  Rk  has  already  been  calculated  need  not  be  redone. 

We  have  written  computer  codes  to  handle  the  above  three  elements 
which  go  into  the  l  calculation.  These  codes  employ  all  the  symmetry 
relations  we  have  discussed.  The  calculation  of  the  matrix  elements  V  in 
Eq.  (3.7)  requires  a  major  effort.  We  have  used  very  efficient  quantum- 
chemistry  computer  codes  for  calculating  two-electron  integrals  over  Gaussians. 
By  fitting  each  local  orbital  to  contracted  sets  of  Gaussians,  we  are  able 
to  write  the  matrix  elements  V  in  terms  of  the  Gaussian  integrals. 
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Tlu'  transformation  in  Eq.  (3.7)  is  then  carried  out  using  quantum  chemistry 

integral  transformation  programs  for  configuration  interaction  calculations. 

*  -> 

We  use  only  up  to  next  nearest  neighbors  in  the  sum  over  R  ,1^  -*-n 

Kq.  (3.3)  for  our  TiO  calculation  since  we  find  adding  more  shells  of  atoms 

do  not  affect  the  major  results  of  the  correlation  calculation.  We  have 

onlv  gotten  results  for  TiO  at  k=0  so  far,  but  these  indicate  that  the  I'  (d) 

_(p)  points  move  away  from  one  another  as  Kunz  expected.  We  will  need  to 

calculate  other  points  in  the  zone,  however , before  we  can  make  out  the 

positions  of  the  correlated  bands.  Our  results  indicate  that  the  diagonal 
v2  >  -> 

elements  of  ]  (k,w)  give  a  qualitative  picture  of  the  correlation  effects 
even  without  diagonalizing  Eq.  (3.1). 

After  we  complete  the  TiO  studies,  we  will  use  the  procedures  we  have 
developed  to  study  correlation  in  NiO  and  FeO.  For  these  studies  we  will 
need  self-consistent  energy  bands  and  wave  functions  in  a  Gaussian  basis. 

We  have  implemented  bands  structure  codes  here  due  to  Wang  and  Calloway'7  and 
have  modified  them  to  accept  two  atoms  per  unit  cell  for  these  calculations. 

B.  Electron  Energy  Loss  Structure  of  Clean  and  Oxidized  Transition 
Metals. 

The  characteristic  energy  loss  structure  of  a  material  provides  a  direct 
measure  of  the  contributions  from  interband  transitions  and  collective  modes 
sucli  as  plasmons  to  the  electronic  excitation  spectrum.  The  loss  structure 
is  altered  by  changes  in  bulk  properties  and  also  by  surface  modification  or 
adsorption,  particularly  certain  features  which  are  quite  sensitive  to  adsorption. 
Ln  additon,  since  the  loss  structure  provides  a  direct  measure  of  plasmon  energies 
and  intensities,  one  can  use  it  to  test  the  basic  understanding  of  these  modes. 
These  tests  are  especially  important  for  systems  with  strong  crystal  potentials, 
such  as  the  transition  metals,  where  the  plasmon  behavior  is  not  expected  to 


follow  the  simple  f ree-electron  result. 
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Unfortunatelv  there  i  ••  y^induve  disagreement  between  the  experimental 
result  Liiemselves,  and  in  addition  the  results  have  been  given  different 
interpretations  so  that  very  different  energies  for  plasmons  and  interband 
transitions  have  been  reported.  To  clean  up  the  confusion,  we  have  under¬ 
taken  a  caret ul  reexamination  of  the  loss  structure  for  certain  transition 
metals  by  analyzing  a  new  series  of  measurements  made  by  Keith  Legg  in  our 
department.  These  measurements  are  for  clean  and  oxidized  polycrystalline 
iron  samples  and  also  clean  and  oxidized  polycrystalline  titanium  samples  for 
comparison.  This  analysis  fits  into  the  objectives  of  this  program  since  it 
has  allowed  us  to  clearly  identify  plasmon  structures  and  has  provided  data 
for  comparison  with  the  more  basic  theoretical  calculations  we  have  undertaken. 

The  characteristic  loss  structure  measurements  fall  Into  two  classes: 

measurements  of  losses  suffered  by  electrons  with  energies  of  many  kilovolts 

as  they  are  transmitted  through  thin  films  and  measurements  of  losses  suffered 

by  electrons  with  energies  around  lOOeV  as  they  are  reflected  from  the  material. 

The  former  experiments  are  dominated  by  electrons  which  do  not  undergo  elastic 

scattering  while  the  latter  are  more  highly  resolved  but  contain  electrons 

which  have  been  elastically  scattered  at  least  once.  The  two  types  of  spectra 

9  10 

for  evidently  clean  Ti  surfaces  are  shown  in  Fig.  1  ’  .  There  are  large 

qualitative  differences  between  the  spectra  at  energy  losses  below  20eV, 
especially  between  the  lower  reflection  spectrum  and  the  transmission  spectrum. 
The  data  from  our  measurement  for  clean  Ti,  which  is  a  reflection  measurement, 
i.s  shown  in  Fig.  2a.  It  is  more  highly  resolved  than  the  reflection  data  in 
Fig.  la,  but  it  is  qualitatively  in  agreement  with  it. 

Our  analysis  of  the  clean  Ti  data  is  as  follows.  The  peak  around  5eV  is 
undoubtedly  a  surface  peak  since  it  is  quite  sensitive  to  surface  contamination 
and  since  it  is  less  pronounced  for  higher  primary  electron  energies  as  shown 


in  fig.  la.  The  latter  behavior  is  consistent  with  a  surface  excitation  since 


Fig.  lb  Transmission  energy 
loss  spectrum  (full 
curve)  for  a  titanium 
film. 
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the  mean  free  paths  of  the  primary  electrons  are  greater  at  higher  energies 

which  means  the  bulk  signal  would  be  enhanced  with  respect  to  the  surface. 

Wo  have  identified  this  peak  tentatively  as  a  surface  plasmon  in  agreement  with 

9 

Simmons  and  Scheibner  (SS)  since  collective  modes  should  be  more  prominent 
at  this  energy  relative  to  interband  excitations.  However  we  need  to  confirm 
this  identification  with  further  measurements  and  analysis  of  how  the  intensity 
of  this  peak  changes  as  the  scattered  electrons  deviate  from  the  specular 
direction.  An  intensity  maximum  at  an  angle  away  from  the  specular  direction 
would  indicate  that  the  peak  is  a  surface  plasmon. ^  The  prominent  peak  at 
about  lOeV  is  the  volume  plasmon.  The  structure  extending  below  this  peak  for 
about  20eV  probably  contains  the  second  volume  and  surface-volume  contributions 
postulated  by  SS,  though  our  higher  resolution  data  shows  this  structure  to  be 
less  prominent  than  theirs.  However,  we  also  find  a  double  peaks  at  about  35eV 
which  also  appears  in  the  appropriate  place  in  the  iron  data,  and  we  identify  this 
as  the  split  (  ,  -^)  Ti  3p  interband  excitation. 

The  loss  spectrum  changes  greatly  upon  adsorption  of  only  a  monolayer  of 
oxygen.  Most  significantly,  the  surface  peak  disappears  completely  and  new 
losses  appear  around  20eV  and  45eV.  The  sensitivity  of  the  spectrum  to  oxygen 

points  out  one  of  the  causes  for  differences  between  the  measurements  of  others 

—8 

workers.  Only  carefully  cleaned  Ti  specimens  in  a  vacuum  of  10  Torr  or  less 
exhibited  the  true  clean  surface  loss  spectrum.  We  identify  enhancement  in 
loss  structure  at  about  20eV  as  due  to  the  0  2s  interband  excitation.  A  similar 
structure  appears  in  the  same  place  in  oxidized  iron.  The  peak  at  about  45eV 
is  more  mysterious.  Transmission  measurements  for  Ti  (presumably  partly 
oxidized)^  and  TiO^  show  this  structure  quite  prominently  in  both  the  loss 
spectrum  and  the  imaginary  part  of  the  dielectric  function  derived  from  a 
Kramer s-Kronig  analysis.  Its  appearance  in  means  it  is  not  a  collective  mode 
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and  must  be  an  interband  transition.  A  similar  peak  does  not  appear  in 

the  oxidized  iron  spectrum  (see  Fig.  3b),  and  there  are  no  interband  losses 

expected  in  this  region.  No  strong  explanation  has  ever  been  offered  for 

ibis  peak,  though  it  could  be  a  double  interband  excitation  or  an  excitation 

to  a  higher  conduction  band.  It  seems  possible  that  the  structure  is  connected 

to  the  metal  3p  peak  since  it  appears  about  9eV  above  this  peak  for  both 

13 

titanium  and  vanadium.  It  could  be  a  double  excitation  involving  the 
metal  3p  and  0  2p  electrons,  but  it  would  be  surprising  that  such  a  process  would 
be  stronger  than  the  3p  excitation  alone.  We  certainly  need  further  theoretical 
analysis  here  and  have  undertaken  an  investigation  of  the  calculated  loss 
spectrum  in  this  region. 

But  what  about  the  differences  between  the  transmission  spectrum  for  Ti  in 

Fig.  lb  and  our  reflection  spectra  in  Fig.  2?  Fig.  lb  most  resembles  Fig.  2b 

because  the  transmission  sample  undoubtedly  had  an  oxide  layer.  However  this 

does  not  explain  the  difference  in  energy  between  the  volume  plasmon  peaks  in 

the  two  experiments,  and  this  difference  is  important  in  an  analysis  of 

collective  modes  in  transition  metals  since  the  transmission  spectrum  places 

the  plasmon  near  its  free  electron  value  while  the  reflection  spectrum  does  not. 

Fig.  la  and  also  reflection  experiments  on  vanadium  done  at  several  primary 
13 

electron  energies  indicate  that  the  plasmon  peak  moves  to  higher  energy  as  the 
primary  energy  increases.  These  results  imply  that  the  differences  between  the 
plasmon  peaks  in  the  two  spectra  may  be  due  to  the  great  difference  in  primary 
energies  used.  A  primary  energy  increase  would  lessen  the  lateral  spread  of 
the  electron  beam  due  to  multiple  scattering  for  a  given  target  thickness. 
However,  in  the  reflection  experiment  as  the  primary  energy  increases,  the 
depth  of  the  target  probed  increases  also.  We  can  do  a  crude  calculation  based 
on  the  plural  scattering  theory  which  should  give  us  an  estimate  of  the  lateral 
spread  for  electrons  transmitted  through  a  300  8  film  at  33  keV  and  for  reflected 
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elections  at  lOOeV  which  no  through  about  20  X.  The  calculation  shows  that 
the  transmitted  electrons  have  a  spread  in  momentum  transfer  of  several  times 
kj .  Thus  even  though  the  spectrum  in  Kig.  lb  bad  an  analyzer  acceptance  angle 
tor  scattered  electrons  of  only  about  5  mrad.^  which  corresponds  to  a  momentum 
transler  of  less  than  .3  i;  for  Ti,  there  must  be  a  groat  deal  ol  scattering  u! 
electrons  with  larger  momentum  transfers  back  into  this  region.  The  dispersion 
ot  the  plasmon  is 


U),  (q)  =  +  iOfi/m)q^  (3. 10) 

1  P 

in  an  RFA-like  theory.  Since  u  is  positive,  dispersion  increases  the  plasmon 
energy  and  would  result  in  the  effect  noted  for  the  transmission  experiment  as 

14 

well  as  significant  peak  broadening.  For  A>.  the  dispersion  is  as  much  as  lieV, 
so  it  is  possible  to  have  an  average  peak  shift  of  several  eV.  We  thus  conclude 
that  the  peak  shift  in  the  transmission  experiment  is  due  to  dispersion  and  that, 
in  this  case,  the  reflection  experiment  peak  is  closer  to  the  true  volume  plasmon 
energy. 

Uur  loss  spectra  for  c  ean  and  oxidized  polycrystalline  iron,  shown  in 
Fig.  3,  have  some  similarities  with  the  titanium  spectra.  The  volume  plasmon 
peak  at  15eV  is  well  below  the  free  electron  value  of  about  30eV,  and  the  spin- 
orbit  split  Fe(3p)  transition  is  clearly  visible  in  the  appropriate  place  at 
about  55eV.  However,  in  contrast  to  titanium,  the  interband  excitations  are 
much  stronger  in  the  iron  spectra.  The  region  below  lOeV  is  dominated  by 
interband  excitations  from  the  iron  d-bands  which  completely  cover  the  effect 
of  the  surface  plasmon.  As  the  oxide  grows  a  pronounced  0  2s  peak  appears  at 
about  19eV  and  grows  larger  than  the  volume  plasmon.  The  structure  below  lOeV 
is  also  changed  due  to  the  growth  of  the  0  2p  peak  and  some  charge  transfer 
from  the  surface  d  bands.  The  peak  at  about  45-50eV  does  not  appear  in  the 
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ox  id i zed  loss  spectrum  as  it  does  ior  Ti  and  V  which  is  consistent  with  the 
idea  t hat  it  is  connected  with  the  metal  3p  band.  However,  it  does  not  appear 
beyond  the  Ke  ip  peak  either  in  Fig.  3b.  The  inelastic  electron  transmission 
loss  data  tor  iron^  differs  from  the  reflection  data  in  a  similar  manner  Lo 
that  from  titanium.  Here  again  we  believe  the  multiple  scattering  caused  as 
the  electrons  traverse  the  sample  mixes  in  plasmons  from  a  large  range  of 
momentum  transfers  and  results  in  a  shift  of  the  plasmon  peak  to  higher  energy. 

We  are  presently  calculating  the  band  structures  for  both  iron  and  iron  oxide 
with  the  intent  of  making  a  comparative  study  of  how  band  structure  and  transition 
probabilities  affect  the  loss  processes. 

C .  Effects  of  the  Dynamically  Screened  Coulomb  Interaction  on  Bose 
Condensation  and  Superconductivity. 

We  have  recently  extended  our  investigation  of  the  dynamical  response 
of  electron  in  transition  metals  to  the  conditions  of  Bose  condensation  and 
superconductivity.  The  usual  BCS  theory  of  superconduct ivitv  13  approximates 
the  term  involving  the  coulomb  interaction  as  a  constant,  but  the  inclusion  of 
the  dynamical  interaction  allows  the  study  of  a  rich  variety  of  effects  caused 
hv  plasmons  anil  ox<  itons.  Up  to  recently  ideas  about  the  existence  of  strong 

effects  of  this  type  were  merely  speculative,  but  some  experimental  evidence 

,  r  16,17 

now  suggests  that  the  effects  may  exist  for  some  systems.  Also, 

calculations  indicate  that  these  effects  should  be  enhanced  in  two  dimensional 

systems  such  as  surfaces  or  inter f aces . ^ ^  All  this  work  improves  the 

prospects  of  finding  condensates  with  unusual  properties  or  high  critical 

temperatures. 

We  have  been  exploring  ways  to  add  more  accurate!''  the  band  structure 
et tints  to  the  calculations  of  the  properties  of  these  Bose  Einstein  systems. 

Local  field  effects  are  assumed  to  significantly  affect  the  critical  temperatures 
of  certain  excitonic  systems,  although  there  is  disagreement  about  whether  the 
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lomperature  is  enhanced*"  or  depressed.  These  arguments  are  based  on  a 
qualitative  analysis  and  a  more  quantitative  analysis  for  specific  two-dimensional 
and  three-dimensional  systems,  such  as  we  have  undertaken,  would  be  quite  useful. 
We  are  exploring  the  properties  ot  systems  witli  the  Hamiltonian, 


H  = 


V  i:k  ctk  ck  +  A  y  l  v°  c+  <•+  c.  c  . 

■  .  n  nii  no  2  ,  L  pa  at  k-qi  p+qo 

n,k.  p,k  ai  1  1 

q 


(3.1  0 


where  |>=np,p;q  -  n  ,<| ; etc . . .  and  we  have  suppressed  the  indices  of  the  Coulomb 

•  O  k 

matrix  element  V  .  The  index  n  refers  to  band  n,  E  is  a  band  energy, 

n  P  •  •  •  n  ° 

and  all  matrix  elements  are  over  Bloch  states.  The  anomalous  Green’s  function 
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in  tin*  Cor'kov  ladder  approximation  would  then  be  at  T=T 


F 
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(p,iuip) 


(q,iiiq) 


x  '■  n  (p  +  q  » 


iS2  •+  i:  ;  ) . 

P  q 


(3.11) 


where  mp  =  tiT(2p+l)  and  cuq  =  2nT  q,  with  p  and  q  integers.  The  quantity 

V  (  q  ,  ifi  )  is  the  effective  interactions  between  electrons  in  the  n'th 
nn  '  q 

band,  and  we  have  assumed  that  intraband  interactions  dominate  interband 
interactions.  The  effective  interaction  may  be  written 
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where  the  parentheses  (...)  indicates  the  diagonal  matrix  element  of  s  he 
product  of  Vq  and  >_  *.  If  we  were  to  neglect  local  field  effects,  then 
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and  we  could  use  the  usual  dielectric  function  in  the  analysis  of  Eq .  (3.11). 

.his  analysis  has  been  done  for  the  case  of  weak  -  coupling  superconductors 

wilh  i  vi'iv  simple  plasmon  pole  approximation  for  general  surface  and  bulk 
I  H 

•vsi  mis  and  a  more  complicated  screening  function  for  a  Si(lll)-SiO 
1  9 

interlace  system.  We  plan  to  study  the  ini luence  ol  local  field  effects  on 
tiic  .  ritical  temperature  in  the  weak-coupl ing  case  for  some  transition  metal 
oxides.  Since  the  plasmon  or  exciton  mediated  coupling  may  not  be  small  at 
all,  we  will  also  look  at  the  strong-coupling  case  for  these  systems  using  the 
theory  of  Kllashberg.^ 

ib 

Measurements  bv  Witteborn  and  Fairbank  have  established  the  existence 

•■I  an  anomalously  weak  electric  field  outside  the  surface  of  copper  at  low 

' e",peratures.  Analyses  have  shown  that  one  should  expect  a  combined  electric 

iii‘ld  due  to  gravitationally  induced  strain  fields  and  patch  fields  of  the 

order  19  h  V/m,' *  but,  in  fact,  the  electric  field  is  only  about  10  ^  V/M 

below  4 . 5°K.  Above  4.5°K  a  phase  transition  seems  to  occur, ^  and  the 

electric  field  outside  the  copper  surface  is  in  the  expected  range.  Many 

analyses  which  suggest  mechanical  or  electronic  explanations  of  this 

phenomenon  have  been  put  forth,  but  perhaps  the  most  likely  and  most  intriguing 

explanation  is  that  the  surface  electrons  couple  and  obey  Bose-Einstein 

23 

statistics  at  low  temperatures.  They  then  form  a  Bose  condensate  below 

'i.s'k.  The  coupling  mechanism  would  be  something  other  than  phonons,  perhaps 

►  ■X.  it  oits  or  plasm  ons.  The  possibility  of  superconductivity  below  the 

transition  temperature  is  also  being  explored.  Since  the  copper  surface  in 

O  16 

the  experiments  is  covered  with  an  oxide  of  some  20  this  shielding  effect 

must  taki'  place  in  or  near  the  oxide.  Now  very  recent  experiments  on 
indicate  that  an  exciton  gas  formed  in  this  material  by  laser  excitation  obeys 
Bose-'. inst e i n  statistics.  This  may  be  the  mechanism  we  need  to  explain  the 

1 

1 
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unomaLous  shielding  effect  seen  by  Witteborn  and  Fairbank  since  their  copper 
oxide  was  most  likely  reduced  and  could  have  been 

Madey  and  Hanni  have  made  an  estimate  of  the  density  of  surface  bosons 
needed  to  provide  the  observed  shielding  of  the  ambient  electric  field  above 
4.5°K.  There  result  is  about  lO"*'3  bosons/cm3  which  is  in  the  low  density 

region  where  exciton  or  plasmon  coupled  superconductivity  is  expected  to  be 

,  .  18,19 

enhanced . 


However,  if, for  example,  one  assumes  the  excitons  in  Cu^O  behave  like 


an  ideal  Bose  gas,  the  volume  condensation  temperature  is 


.  17 


2/., 
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g6. 2x10 
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(3.13) 


where  mQ  is  the  f ree-electron  mass,  m  is  the  particle  mass,  g  is  the  degeneracy 

_3 

factor,  and  n  is  in  units  of  cm  .  For  Cu^O,  m/m^  3  and,  since  the  most 

Likely  exciton  state  to  obey  Bose-Einstein  statistics  is  the  triplet,3^  g=3. 

11  3/2 

Assuming  a  volume  density  of  n=(10  )  ,  the  result  from  Eq.  (3.13)  is 

T  o 

c  =  .329  K.  However,  we  are  dealing  with  an  essentially  two-dimensional 

24 

system  for  which  Ginzburg  and  Kirshnits  estimate  the  transition  temperature 


for  an  ideal  Bose  gas  should  be 


T 

c 


dn1/3/ln(Adn)  T 


C00 


(3.14) 


where  d  is  the  film  thickness  and  A  is  the  film  area.  With  d=20  X  and 
A  100  cm"3,  =  2.32x10  3  Tc  =  .76x10  3  °K.  This  temperature  is 

CD 

considerably  below  the  observed  transition  temperature.  As  Ginzburg  and 
24 

kirshnits  observe,  Eq.  (3.14)  may  underestimate  the  transition  temperature 
since  superconducting  films  of  thickness  d  ^  30  X  and  T^  =  T^  exist. 
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Unf ortiinately  the  ideal  gas  result,  Kq.  (3.13), is  probably  an  upper  limit  to 
T  .  For  finite  exciton  coupling  constant  A  the  critical  temperature  would 
be  lower  and  still  well  below  the  transition  temperature  observed.  This 
discrepancy  is  a  serious  defect  in  explanations  of  the  anomalously  weak 
electric  field  in  terms  of  screening  by  a  Bose  condensate. 

23 

A  review  of  the  screening  calculations  of  Madey  and  Hanni  indicates  that 

assuming  an  ideal  Bose  gas  composed  of  coupled  free  electrons  may  overestimate 

the  screening.  The  electric  field  in  the  surface  region  due  to  the  system  of 

23 

substrate  and  surface  electrons  and  ions  is 


1  _9p^  ma 

e  3z  e 


(3.14) 


where  m  is  the  chemical  potential  of  the  electrons  and  a  is  the  acceleration 
due  to  gravity.  Now  since  p  varies  mainly  through  the  change  in  the  surface 
density  ng  , 


,  3n 

£p  _  a_p  _ s 

3z  3n  3z 
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4Tie^d 


(3.15) 


The  last  expression  involves  an  estimate  of  the  variation  of  n^  due  to  the 

substrate  electric  field  which  involves  both  patch  fields  and  gravitationally- 

induced  lattice  compression  fields.  The  quantity  - ^ —  (-^— )  is  then  the 

hue  d  s 

shielding  factor.  Now,  for  an  uncondensed  ideal  Bose  gas  in  two-dimensions, 
u  is  defined  by  the  expression 

n  -2  fd2Jo3(E4V/2V,0  i 

%  ~  2^  J  d  k\e  °  0  -1 

where  B  =  1/k  T  and  E  is  the  ground  state  energy.  The  integral  in  Eq.  (3.16) 

D  O 


(3.16) 
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is  analytic  in  the  two-dimensional  case  and  the  result  is 
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-1 


In 


1  -  exp(-2iT'fi“  n  6/m  ) 
s  o 
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(3.17) 


ihe  electric  field  at  the  surface  state  is  then 
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It  is  (it).  (3.18)  that  Madey  and  Hanni  use  to  fit  the  experimental  data  of 

_6  pp  _2 

Witteborn  and  Fairbank  with  an  c  ^  10  V/m  and  n  'v  10  cm  . 

z  s 

We  can  include  the  band  structure  effects  on  the  excitonic  gas  in  Cu^O 
as  before ,  witli  an  effective  mass  replacing  the  electron  mass  in  Eq.  (3.16). 
Ihe  degeneracy  factor  g  should  also  multiply  the  integral  in  Eq.  (3.16),  and 
the  result  for  u  would  be 


=  if1  In 


1  -  exp(-2Trh 


ngB/gm) 


+  E 


(3. 16’) 


Since  g=3  and  m  'v  3mo  for  the  excitons  in  C^O,  we  can  now  fit  the  experimental 
data  with  n  '>  10^  cm  This  value  of  n  gives  T  =  3.3°K.  The  value 

s  s  Co<i 

r  -  10  V/m  for  the  substrate  fields  may  also  be  low  due  to  a  larger  patch 
field . 

It  is  obvious  from  these  estimates  that  we  must  make  a  more  careful 
examination  of  the  shielding  due  to  the  exciton  gas  in  C^O  by  including  in 
our  calculations  band  structure  effects,  exciton  interactions,  and  the  effects 
of  electron-hole  coupling.  The  condensation  temperature  for  an  interacting  gas 
of  excitons  should  also  be  studied.  We  will  also  explore  the  possibility 
that  the  condensed  phase  is  superconducting.  Finally,  the  excitonic  gas  in 
Cu^O  was  obtained  by  laser  excitation  in  the  experiments;  we  must  study  the 
possibility  that  a  surface  layer  of  excess  electrons  can  also  form  the  gas. 


I’li hlications  and  Presentations  • 

The  tol lowing  papers  are  being  prepared  for  publ icat ion: 

1.  "The  Effects  of  Correlation  on  the  Band  Structure  of  TiO". 

2.  "Analysis  of  Knergy  Loss  Spectra  of  Clean  and  Oxidized  Titanium  and 
Iron,"  with  K.  0.  Legg. 

The  following  paper  is  scheduled  for  presentation: 

1.  "Characteristic  Energy  Loss  Spectra  of  Clean  and  Oxidized  Titanium 
Surfaces,"  with  K.  0.  Legg;  March  Meeting  of  the  APS,  Phoeniz,  Az . 
(198i).  An  abstract  for  this  paper  is  attached  to  this  section. 
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session  in  which  paper 
should  be  placed: 

Energy  Loss  Spectroscopy 


Characteristic  Energy  Loss  Spectra  of  Clean  and 
Oxidized  Titanium  Surfaces,  M.W.  Ribarsky  and  K.O.  Legg, 
Georgia  Institute  of  Technology.* — We  present  recently 
taken  energy  loss  spectra  for  clean  and  oxidized 
titanium  and  analyze  the  loss  structures.  The  use  of 
different  primary  electron  energies  and  electron  exit 
angles  and  comparisons  with  optical  absorption  data  help 
identify  surface  and  bulk  plasmons  and  interband  tran¬ 
sitions.  Our  interpretation  of  the  spectra,  especially 
for  the  oxide,  differs  from  earlier  interpretations,  and 
we  discuss  the  implications  of  this  for  our  view  of 
characteristic  losses  in  titanium  oxide  and  for  other 
processes  occuring  at  the  surface.  We  also  analyze  the 
plasmon  structure  in  the  light  of  various  theoretical 
models. 
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V.  Personnel 

The  to] lowing  personnel  have  been  partially  compensated  by  funds 
provided  in  this  grant  for  conducting  the  research  during  the  period 
October  1,  1979  to  September  30,  1980. 

Dr.  Martin  W.  Ribarsky  -  Research  Scientist  and  Principal  Investigator 
Mr.  David  Luedtke  -  Ph.D  candidate. 
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V I .  Interactions 

We  had  the  following  interactions  relative  to  this  research  during  the 
period  October  1,  1979  to  September  30,  1980: 

1.  With  Professor  Joseph  Callaway  and  his  group  at  Louisiana  State 
University  about  our  modification  and  use  of  LCAO  self-consistent 
band  structure  codes  developed  by  his  group. 

2.  With  Peter  Feibelman  and  Dwight  Jennison  at  Sandia  Laboratories 
about  calculations  on  transition  metal  oxides  to  be  used  as  a 
basis  for  calculating  electronic  excitation  structure  and  plasmon 
structure. 

3.  With  Dr.  Brett  Dunlap  of  N.R.L.  about  developing  and  using  LCAO-Xa 
cluster  codes  for  the  calculation  of  the  electronic  structure  of 


surface  oxides. 
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